Electrode Polarization Effects in Broadband Dielectric Spectroscopy 
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In the present work, we provide broadband dielectric spectra showing strong electrode polarization 
effects for various materials, belonging to very different material classes. This includes both ionic 
and electronic conductors as, e.g., salt solutions, ionic liquids, human blood, and colossal-dielectric- 
constant materials. These data are intended to provide a broad data base enabling a critical test of 
the validity of phenomenological and microscopic models for electrode polarization. In the present 
work, the results are analyzed using a simple phenomenological equivalent- circuit description, in- 
volving a distributed parallel RC circuit element for the modeling of the weakly conducting regions 
close to the electrodes. Excellent fits of the experimental data are achieved in this way, demonstrat- 
ing the universal applicability of this approach. In the investigated ionically conducting materials, 
we find the universal appearance of a second dispersion region due to electrode polarization, which 
is only revealed if measuring down to sufficiently low frequencies. This indicates the presence of a 
second charge-transport process in ionic conductors with blocking electrodes. 
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1. INTRODUCTION 

Electrode polarization (EP) leading to blocking elec- 
trodes (BEs) is a frequently encountered phenomenon 
when investigating ionic conductors with impedance or 
dielectric spectroscopy. This effect, giving rise to giant 
values of the dielectric constant and a strong drop of con- 
ductivity towards low frequencies, arises when the ions 
arrive at the metallic electrodes and accumulate in thin 
layers immediately beneath the sample surface forming 
a so-called space-charge region. It occurs in solid-state 
electrolytes as well as in aqueous solutions, ionic liquids, 
and in many biological systems (e.g., [H4l8|). EP in ionic 
conductors often hampers the determination of the in- 
trinsic dielectric properties but it also is of high technical 
relevance, e.g., in the design of fuel cells or double layer 
capacitors. 

EP also plays an important role in the dielectric spec- 
troscopy of electronic conductors (e.g., [l9l - [3ll |). In semi- 
conducting samples it is caused, e.g., by the formation 
of Schottky diodes at the electrode-semiconductor inter- 
face. In this case the blocking of the electrode arises from 
the depletion layer at the semiconductor/metal contact 
in the blocked state of the diode. Recently these effects 
came into the focus of interest due to their presumed 
importance in the so-called " colossal dielectric constant" 
(CDC) materials [U H H H, ■ 

The many suggestions made to account for electrode ef- 
fects in dielectric measurements can be classified into two 
categories, namely, (i.) avoiding the electrode polariza- 
tion by specialized measurement techniques or (ii.) mod- 
eling these contributions, e.g., by equivalent circuits. The 
first approach, using, for example, four-electrode meth- 
ods or varying the electrode distances (e.g., fll. l33l. l34j ) . 
quite generally implies considerable additional experi- 
mental efforts compared to standard dielectric measure- 
ments. Thus, often the second approach is adopted and 
numerous, partly quite sophisticated phenomenological 
and microscopic models have been proposed to model 
EP effects in dielectric spectra (e.g., [IMl H3, IH- 



|43|). In dielectric studies, one often deals with situations 
where the detected EP effects are not in the focus of in- 
terest but still cannot be ignored since they have to be 
taken into account for an unequivocal determination of 
the intrinsic properties. In such cases, for a straightfor- 
ward modeling of EP effects, simple equivalent circuits 
are usually employed, which are assumed to be connected 
in series to the intrinsic bulk contribution of the sample 
material. The most common ones are a parallel RC cir- 
cuit [l(| Oil (22|, [23, [44 . rial o r a so-called constant phase 
element (CPE) [J 07 0G3 M, HI 0, the latter be- 
ing mainly applied for modeling ionic conductors. Both 
models work relatively well if only the onset of EP ef- 
fects is seen, arising close to the low-frequency boundary 
of the covered frequency range. However, in our experi- 
ence both approaches have certain drawbacks when ana- 
lyzing broadband dielectric spectra, where contributions 
from EP are observed over several frequency decades. In- 
stead we found that a distributed RC equivalent circuit 
(referred to as DRC circuit henceforth) , based on a Cole- 
Cole (CC) distribution of relaxation times [2l|, HtHMJ , is 
best suited to describe EP. This is true for such different 
materials as aqueous solutions, biological systems, solid- 
state electrolytes, ionic liquids and melts, and all kinds 
of electronic conductors. 

The purpose of the present work is twofold: i.) To 
provide broadband experimental spectra with strong EP 
effects in a variety of very different materials that may 
be helpful for future tests of models of EP (for this pur- 
pose, the experimental data are provided in electronic 
form in the supplementary information, ii.) To demon- 
strate the universal applicability of the DRC circuit for 
the modeling of EP effects. 



2. MATERIALS AND METHODS 

Low-frequency measurements up to v ~ 3 MHz were 
performed using frequency-response analyzers (Novocon- 
trol Alpha-A and Schlumberger 1260 combined with the 
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Chelsea dielectrics interface) and an autobalance bridge 
(Hewlett-Packard HP4284A). Here the sample always 
was placed in a parallel-plate capacitor. Data in the 
higher frequency ranges were collected using a coax- 
ial reflection technique, employing the Agilent Impe- 
dance/Material Analyzer E4991A (1 MHz - 3 GHz) and 
the autobalance bridge 4294 A (40 Hz - f 10 MHz). Partly, 
additional measurements beyond 3 GHz were performed 
via an open-end coaxial reflection technique and a coax- 
ial transmission technique using the network analyzers 
Hewlett-Packard 8510 and Agilent E8363B. Further de- 
tails on the experimental techniques can be found in refs. 
d, [H, HH, HH . Cooling and heating of the samples were 
achieved by a N2-gas cryostat (Novocontrol Quatro) and 
various home-made ovens. 

In the present publication we provide broadband di- 
electric data for the following samples (all concentrations 
in mol percent): 

i. ) LiBr/water and LiBr/glycerol solutions with the 
concentrations c — 0.018 % and c = 0.73 %, respectively. 
For the preparation of the solutions, 54 wt% aqueous 
LiBr solution (Sigma Aldrich) , Glycerol (AnalaR Norma- 
pur, 99.5%), H 2 (Merck "Ultrapur"), and LiBr powder 
(Sigma Aldrich, purity 99.999%) were used. 

ii. ) LiCl/water solution with c = 17 % (Sigma Aldrich, 
purity > 99%). 

iii. ) Human whole blood (hematocrit value of 0.39) 

m. 

iv. ) An aqueous protein solution: lysozyme with c = 
0.009% (5 mmol/1) (lysozyme: Sigma Aldrich (Fluka)). 

v. ) The glass- forming ionic melt 
[Ca(N0 3 )2]o.4[KN0 3 ]o.6 (CKN) (see ref. Q for 
preparation details). 

vi. ) The ionic liquid l-butyl-3-methylimidazolium- 
bis(trifiuoromethylsulfonyl)imide (BMIM-BMSF) (Io- 
litec, purity 99%) 

vii. ) Single crystalline La2Cu04, a parent compound 
of high- T c superconductors (see refs. 0,[2l| for details). 

viii. ) Single crystalline La 15 / 8 Sr 1 / 8 Ni04, a charge- 
ordered CDC material 

ix. ) Polycrystalline Pr 2 /3Cu3Ti40i2, a CDC material 



EQUIVALENT-CIRCUIT MODELS FOR 
ELECTRODE POLARIZATION 



The CPE is often employed to account for EP in ionic 
conductors [f| 0, E3, 0, US Ef3|. Its admittance is de- 
fined by Y* = A(iv) n with an exponent n < 1 and the 
prefactor A. A theoretical rationale for the CPE can be 
obtained when considering a fractional geometry of the 
electrode/sample interface [4j, [5J. A more intuitive mod- 
eling of EP effects is provided by a parallel connection of 
a resistance Re and a capacitance CV, connected in series 
to the bulk impedance Z B |3 H El H, IH ■ Here 
Ce models the high capacitance of the thin blocking lay- 
ers and i?E accounts for their high resistance (Re = oo 



for complete blocking). This equivalent circuit leads to a 
total impedance of 



Ztotal — Ze 



Ry, 
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with te = ReCe and u> = 2ttv the circular frequency. 
Assuming, e.g., a frequency-independent bulk resistance 
and capacitance (i.e., another RC circuit) for the bulk 
response, this results in so-called Maxwell- Wagner re- 
laxation. It leads to a frequency dependence of the to- 
tal complex capacitance (and thus of the permittivity 
e* = e' — i e") that is fully equivalent to the Debye relax- 
ation law [54j descri bing the relaxational response of an 
ideal dipolar system [32[ [55[ : 
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r denotes the relaxation time (not identical with te |32[) 
and Ae = e s — the dielectric strength. e s and 
are the limiting values of the real part of the permittiv- 
ity for frequencies well below and above the relaxation 
frequency v re i ax = 1/(2ttt), respectively. The last term 
in eq. @ was included to account for the contribution 
of dc charge transport to the dielectric loss e" cx a'/v, 
with Ode the dc conductivity. Typical spectra of e' and 
the conductivity a' arising from the Maxwell- Wagner re- 
laxation mechanism, described above, are schematically 
shown by the solid lines in Figs. QJa) and (b), respec- 
tively. It should be noted that in case of the validity 
of eq. (Q} , an intrinsic dipolar relaxation and a Maxwell- 
Wagner relaxation caused by EP cannot be distinguished 
based on the measured spectra alone. 

However, we found that in many cases the above ap- 
proach does not lead to a satisfactory description of the 
experimental data, especially if the spectra extend to fre- 
quencies sufficiently low to reveal more than just the on- 
set of the EP effects. Then the spectral features are of- 
ten found to be more smeared out than expected for this 
simple equivalent circuit. Interestingly, intrinsic relax- 
ations, e.g., in dipolar glass forming liquids [5l|, usually 
exhibit a similar behavior, i.e. a broadening of the ob- 
served spectral features. In these cases, the use of an 
empirical extension of the Debye function, the Havriliak- 
Negami function (56|, is well established to account for 
the observed deviations from Debye behavior: 
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The above function with < a < 1 and < (3 < 1 in- 
cludes the special cases with j3 = 1 and a = 0, termed 
Cole-Cole [5]} and Cole-Davidson [HI, [5!| function, re- 
spectively. These two additional parameters lead to a 
broadening of the spectral features, which is commonly 
ascribed to a distribution of relaxation times. While 
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FIG. 1: (a) Dielectric constant and (b) real part of the conduc- 
tivity calculated according to the "Cole-Cole"-case (/3 = 1) of 
eq. for various values of the distribution parameter a. 
a = corresponds to the undistributed case, eq. {]]). For 
the bulk, frequency independent e' and a' were assumed, i.e. 
a simple parallel RC circuit connected in series to the impe- 
dance of eq. U Frames (c) and (d) show corresponding curves 
for a CPE with n = 1 and n — 0.8 instead of a DRC element. 



the Havriliak-Negami and the Cole-Davidson functions 
are purely empirical, the Cole-Cole distribution of r can 
be approximated by the microscopic model of Gauss- 
distributed energy barriers [60]. This formula leads to 
a symmetric broadening (compared to the Debye case) 
of the relaxation peak in e" and, e.g., is often used for 
the description of secondary relaxations of glassy matter 
[U [Hj]- However, we found that a direct application of 
eq. (J3J to Maxwell- Wagner relaxations arising from EP 
cannot account for the experimental data. An alterna- 
tive approach is an extension of eq. ([1]) analogous to eq. 
@, arriving at 



-'total 



Re 
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(4) 



with < a < 1, < j8 < 1. For the Cole-Cole case 
(8 = 1), this equation was already considered earlier (47l — 
150 ] . A generalized Havriliak-Negami version as in eq. (U) 
is also included in the "LEVM" fit programm by Prof. 
J. Ross Macdonald [631 ]. However, until now this DRC 
equivalent circuit has only rarely been used for the mod- 
eling of EP effects [U 13 |64| • As shown in the further 
course of the present work, it seems that the CC version 



of eq. ((4]) provides good fits of EP-dominated spectra in 
a large variety of very different materials. 

It has to be stressed that, whereas a simple RC equiv- 
alent circuit (eq. (J])), connected in series to the bulk, 
produces permittivity curves identical to a Debye func- 
tion (eq. ©), a DRC equivalent circuit based on eq. ^ 
in general produces permittivity spectra that are quali- 
tatively different from those of a Havriliak-Negami func- 
tion, eq. ([3]). For the latter, the relaxation time r in Eq. 
((3]) is assumed to be distributed. In the equivalent-circuit 
case, the corresponding quantity is r = RbCe (Rb de- 
notes the bulk resistance), which determines, e.g., the 
loss pe ak p osition via v v = 1/(2ttt), just like r in eq. 
© [24, 32]. However, the quantity that is distributed in 
the equivalent-circuit case is te = ReCe, thus leading to 
different curve shapes. 

Finally it should be mentioned that the DRC circuit 
for the Cole-Cole case and a parallel connection of a CPE 
and a resistor can produce identical frequency depen- 
dences of the dielectric quantities The DRC circuit 
seems to be the more physical one because fitting the 
experimental data with this model provides direct ac- 
cess to the resistances and capacitances of the blocking 
surface layers. However, as will become clear in the fur- 
ther course of this work, the extended CPE element may 
find an explanation in a seconchslower hopping transport 
process within the BE layers fla ]. 

Figure Q] shows typical spectra of e' (a) and a' (b) cal- 
culated by assuming a DRC equivalent circuit to account 
for the EP. For comparison, spectra using a plain RC 
circuit (solid lines in (a) and (b)) and a CPE (frames 
(c) and (d)), instead of the DRC, are provided. In all 
cases for the bulk response a simple parallel RC circuit 
was assumed, corresponding to a frequency independent 
bulk e' and a'. For the distribution, the CC case (i.e., 
/3 = 1 and a ^ in eq. ((4]) was chosen. The curves 
for the undistributed RC circuit (solid lines), which are 
identical to Debye behavior, eq. ©, show the typical 
dispersion steps in e' and a' . With increasing a, these 
steps become smeared out for both quantities. In addi- 
tion, the low- frequency plateau of e'{v) develops into a 
weak power law. The low and high-frequency plateaus 
of <j'{v) remain unaffected (however, we found that for 
very high values of a the o'(y) curve can cross those lim- 
its). Comparing this behavior with that of the CPE case 
(Fig. [IJc) and (d)), the most obvious difference is the 
missing low- frequency plateau in a'(y) for the latter (d). 
This would correspond to completely blocking electrodes. 
However, adding a parallel resistance to the CPE leads 
to such a plateau and, as mentioned above, in this case 
both circuits produce identical curves. 

In the following, several experimental spectra being 
governed by strong EP effects, obtained on partly very 
different materials, will be fitted by the models intro- 
duced above. All fits were simultaneously performed for 
e'{v) and o'(y). For the experimental data, statistical 
errors were assumed. 
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RESULTS AND DISCUSSION 
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FIG. 2: (a) Dielectric constant and (b) real part of the con- 
ductivity of CKN at 379 K as function of frequency. The 
lines represent fits using different ways to account for the EP 
effects: equivalent circuits with a CPE (dash-dotted lines), a 
parallel RC circuit, eq. ([T]) (dotted lines), or a DRC circuit, 
eq. Q with /3 = 1 (solid lines), connected in series to the 
bulk. For the latter, a power-law increase in a'{v) and the 
corresponding contribution in e'(y) was used [HEI ]. For com- 
parison, the dashed lines show fit curves where the EP effect 
was treated like an intrinsic relaxation and fitted by a Cole- 
Cole function plus dc-conductivity contribution (eq. ([3]) with 
= 1). 

Figure [2] shows broadband spectra of e'{y) and <x'(z/) 
of CKN at 379 K, extending over 12 frequency decades 
[1, |65j]. In this supercooled ionic melt, EP effects lead to 
a strong increase of the dielectric constant and a strong 
decrease of the conductivity at low frequencies. The in- 
crease of <j'{v) at high frequencies is of intrinsic nature 
and corresponds to an onset of a high-frequency min- 
imum in the dielectric loss, e" oc a' j v. This typical 
feature of glass formers is treated in detail in [65|, |66| . 
where data extending to even higher frequencies are pro- 
vided. The dotted and dash-dotted lines in Fig. [2] show 
fits with the two standard descriptions of EP, namely a 
parallel RC circuit and a CPE, respectively, both con- 
nected in series to the bulk. For the bulk, in addi- 
tion to the usual parallel RC circuit, we have assumed 
a power-law increase in u' (y) and the corresponding con- 
tribution in e '(is) to phenomenologically account for the 



mentioned high-frequency effects [5a ] . The EP-generated 
low- frequency dispersions in e'{v) and a'{y) of CKN are 
clearly too smeared out to be fitted by a simple RC cir- 
cuit (dotted line). In addition, the absence of a clear 
low- frequency plateau in s'(u) is not taken into account 
by this approach. Applying this RC equivalent-circuit is 
mathematically equivalent to fitting the data with an in- 
trinsic Debye relaxation, eq. @. Thus one may expect 
better fits if using a broadened function like the Cole-Cole 
function, eq. §3§ with /3 = 1. However, as demonstrated 
by the dashed lines in Fig. [51 such fits are still unsat- 
isfactory. This is also the case for the Havriliak-Negami 
and Cole-Davidson variant of eq. ([3]) (not shown). 

While the CPE (dash-dotted lines in Fig. ^ does a 
good job in the description of e'(u ), it fails to account for 
the onset of a low- frequency plateau observed in o'iy) 
at v < 0.1 Hz. This plateau, which is often overlooked 
because of a limited frequency range or low ion mobility, 
obviously is due to an incomplete blocking of the elec- 
trodes. It indicates a second conductivity mechanism 
that starts to prevail at v < 10 _1 Hz when the main 
charge-transport mechanism is blocked. It could be due 
to electronic conductivity or to a second ion species with 
much lower mobility, which is not yet blocked at the low- 
est frequencies investigated 14]. However, it seems un- 
likely that the mobilities of the different ion species (K + , 
Ca 2+ , and NOg in the present case) differ by about three 
orders of magnitude as suggested by a comparison of the 
plateaus at v — ► and v k, 10 4 Hz [67]. Another pos- 
sibility is that the ions located within the space charge 
region, close to the electrode, still have some mobility left. 
In line with the model proposed in ref. (l5j , this mobility 
can be assumed to be much lower than in the bulk, but 
not zero, thus explaining the approach of an additional 
low- frequency plateau in a'{y). Then it seems plausible 
that at extremely low frequencies a second, final blocking 
effect should arise, i.e. a'(y) should decrease again when 
the ions within the space charge region finally reach the 
sample boundary. As shown in the further course of this 
work, this is indeed observed in several of the investigated 
materials. 

As demonstrated by the solid lines in Fig. [5J mod- 
eling the BEs by a DRC circuit yields satisfying fitting 
curves of e'{v) and cr 1 (y) in the whole frequency range, 
including the onset of a low-frequency plateau in the con- 
ductivity. The resistance arising from the second charge- 
transport mechanism is taken into account by the resis- 
tor within the parallel RC equivalent circuit, modeling 
the EP effects. We find a = 0.26 and (3 — 1, indicat- 
ing a moderate, symmetrical distribution of relaxation 
times te (Cole-Cole case). Broadband spectra obtained 
at other temperatures in CKN 0, |65[ also can be well 
modeled in this way. In fact, in an earlier work [6~H a so- 
phisticated model that included a DRC circuit (termed 
" ZC" there) to account for EP, was used to obtain excel- 
lent fits of spectra in CKN at four temperatures between 
342 and 361K. 

Figure [3] provides results on BMIM-BMSF, an ionic liq- 
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FIG. 3: (a) Dielectric constant and (b) real part of the con- 
ductivity of BMIM-BMSF at 252 K as function of frequency. 
The lines represent fits assuming different models for the EP 
effects: the sum of two Cole-Cole functions, eq. ([3]) with f3 — 1 
(dashed lines), two parallel RC circuits, eq. (0, connected in 
series to the bulk (dotted lines), and two DRC circuits, eq. 
with /8 = 1, also connected in series to the bulk impedance 
(solid lines). For the latter two cases, the bulk is modeled by 
a parallel RC circuit. 



uid, i.e. a substance solely consisting of cations and an- 
ions that is liquid at room temperature. This class of ma- 
terials has attracted considerable interest in recent years 
due to their possible technical relevance as environment- 
friendly solvents and as new types of electrolytes used, 
e.g., in batteries (H,^. Again, typical EP effects arc 
observed in e'(v) (a) and a'{v) (b). The upward bend- 
ing of the e'{u) curve in the upper-plateau region of the 
main dispersion step below v w 10 Hz clearly indicates 
the presence of a second dispersion step. In Fig. |3£b), 
after the initial decrease of a'(v) below about 3 x 10 3 Hz, 
one may suspect the approach of a second plateau at 
around 10 Hz with amplitude of ~ 10 -6 fi — 1 cm — which 
is followed by a second decrease below about 1 Hz. As 
discussed above, the latter finding may well indicate the 
final blocking of a second charge-transport mechanism. 
Admittedly, the mentioned features are not well pro- 
nounced and rather smeared out. However, we found 
that fits with a single DRC element are not able to pro- 
vide a good description of the present data. Instead, the 



experimental spectra can be very well fitted by assum- 
ing two DRC elements in series to the bulk (solid lines) 
modeling the blocking of the two transport mechanisms. 
Here the bulk is assumed to have frequency- independent 
e' and a', i.e. it is modeled by another (undistributed) 
parallel RC circuit. From the fits, we obtain a — 0.18 
and (3 = 1 for the DRC element, modeling the main 
dispersion step. For the low-frequency one, a pure RC 
element (a = 0, (3 = 1) was sufficient, which may be 
due to the fact that it is only partly visible in the in- 
vestigated frequency range. As indicated by the dotted 
lines in Fig. [31 a combination of two undistributed RC 
circuits (eq. ([J)) in series to the bulk, cannot explain the 
observed frequency dependence. Using instead the sum 
of two Cole-Cole functions (eq. © with (3 = 1; dashed 
lines), i.e. treating the data as arising from two intrin- 
sic dipolar relaxations with different time constants, also 
does not account for the experimental data. 
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FIG. 4: (a) Dielectric constant and (b) conductivity of a 
0.018 % LiBr/water solution at 295 K as function of frequency. 
The meanings of the lines are the same as in Fig. [3] 

As an example for a classical electrolyte solution, Fig. 
|4] provides the frequency dependence of e' and a 1 of a 
0.018 % LiBr/water solution at 295 K. The real part of 
the permittivity (a) reveals two smeared-out dispersion 
steps with points of inflection at around 10 and 10 5 Hz, 
which are caused by EP effects. The conductivity (b) 
also shows two dispersion features. The high-frequency 
plateau in <r'(v) of about 10~ 3 f2 _1 cm _1 corresponds to 
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the intrinsic ion conductivity of this solution, dominated 
by the main charge-transport process (see discussion of 
the CKN data) . This system exhibits a high conductivity 
and an onset of the EP effects at rather high frequencies 
(about 10 5 Hz in a' (v), compared to 3 x 10 3 Hz in BMIM- 
BMSF at 252 K and 10 2 Hz in CKN at 379 K, cf. Figs. 
[3] and [2]) . Probably for this reason, the appearance of a 
second plateau (at about 10 2 Hz and ~ 10~ 5 fi~ 1 cm _1 ) 
is well pronounced in this electrolyte and, in addition, 
the final blocking is well seen (further decrease of u'(v) 
below ~ 10 Hz). Both effects clearly indicate a second 
conductivity mechanism that becomes blocked at lower 
frequencies, as already suspected in the discussion of the 
CKN and BMIM-BMSF data. 

The lines in Fig. [4] again are fits with different em- 
pirical functions to account for the BEs. As demon- 
strated by the dotted lines, two parallel RC equivalent 
circuits (eq. (JTJ), both connected in series to the bulk 
impedance, obviously cannot account for the experimen- 
tal data. Here again the bulk was assumed to have 
frequency- independent e' and a' . The sum of two Cole- 
Cole functions (eq. ([3]) with j3 = 1) leads to much bet- 
ter, but still not perfect fits (dashed lines). In contrast, 
assuming two DRC elements, eq. (fj|, with f3 — 1, con- 
nected in series to the bulk (solid line), provides nearly 
perfect fits. The values of the distribution parameter a 
are 0.32 (for the low frequency EP effect) and 0.28 (for 
the high frequency EP effect). 

As an example for the temperature-dependent devel- 
opment of EP effects, Fig. [5] shows e'{v) (a) and u'{y) 
(b) of a 0.73% LiBr/glycerol solution for various temper- 
atures between 220 K and 410 K. For the lowest temper- 
ature shown (220 K) only the intrinsic a-relaxation of 
glycerol (cooperative reorientation of the dipolar glyc- 
erol molecules [5l{ | ) and the bulk conductivity arising 
from the ionic charge transport can be seen in this fre- 
quency range. This relaxation is characterized by steps in 
e'iv) and in u'(y) (the latter corresponding to the well- 
known loss peaks in e"(y) 51]). The bulk conductiv- 
ity leads to the low- frequency plateau in <r'(v) of about 
10~ 10 f2 _1 cm _1 . Due to the low ion mobility at this 
temperature, the EP effects lie outside the investigated 
frequency range. With increasing temperature, the mo- 
bility rises and the EP shows up as a strong increase 
in e'{v) at low frequencies that finally develops into a 
stcplike curve. In addition, a decrease in the conductiv- 
ity for decreasing frequency arises above 240 K. At even 
higher temperatures, a second, smeared-out step is re- 
vealed in e'(y), accompanied by a weak corresponding 
feature in o'(y). Similar to the behavior in the aqueous 
LiBr solution (Fig. [4|, this further dispersion indicates a 
second charge-transport mechanism as discussed above. 
As demonstrated by the lines in Fig. [5j perfect fits of 
all these spectra are achieved by assuming a DRC circuit 
(eq. ((U) with j3 = 1 for the description of the observed 
EP effects, connected in series to the bulk response. For 
the latter, the a-relaxation is accounted fo r by the Cole- 
Davidson function, eq. ([3]) with a = |5ll . [70| . The 
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FIG. 5: (a) Dielectric constant and (b) conductivity of a 
0.73 % LiBr/glycerol solution at different temperatures as 
function of frequency. The lines are fits with two DRC cir- 
cuits, eq. Q with ft — 1, in series to the bulk to account for 
the EP. For the bulk, an intrinsic dipolar relaxation, described 
by the Cole-Davidson formula and a frequency-independent 
dc conductivity were used (eq. (0 with a — 0) . 



distribution parameters a of the slower EP effect vary 
between 0.1 and 0.33; for the second one we obtained 
values between 0.1 and 0.29. 

Further examples of dielectric spectra, dominated by 
EP are provided in Fig. HI It shows the dielectric proper- 
ties of a variety of partly very different materials, includ- 
ing a highly concentrated salt solution, two biological 
samples, and three electronic conductors. The squares 
represent the results obtained on a 17% LiCl/water so- 
lution at 290 K. Its frequency spectra reveal two strong 
dispersion steps in e'(v) (a) and o 1 (y) (b). As demon- 
strated by the lines, the spectra are well fitted with two 
Cole-Cole DRC elements in series to the sample impe- 
dance, which is modeled by a simple RC circuit. The 
broadening of the low-frequency dispersion step is char- 
acterized by a = 0.22; the faster one could be modeled 
with a pure RC equivalent circuit (i.e. a = 0, /3 = 1). 

The triangles in Fig. [5] represent the broadband data 
of human blood at 290 K [lq . Blood can be consid- 
ered as an aqueous solution of a variety of different salts 
with the addition of a considerable fraction of suspended 
"particles" ■ mainly red blood cells. As treated in detail 
in ref. [18|, at v > 10 5 Hz, those data are dominated 
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FIG. 6: (a) Dielectric constant and (b) real part of the con- 
ductivity of different samples that are affected by EP effects: 
17% LiCl/water solution at 290 K (squares), human whole 
blood at 290 K (triangles), 0.009% Lysozym/water solution 
(diamonds), single-crystalline La2Cu04 (down triangles), sin- 
gle crystalline Lais/sSri/gNK^ (crosses), and poly crystalline 
P r 2/3Cu3Ti4 0i2 (stars). The lines are fits with one or two 
DRC circuits, connected in series to the bulk, which is mod- 
eled by a parallel RC circuit. 



by the well-known f3- and 7-dispersions of blood [TH, [7l[ . 
The /3-dispersion (at about 10 6 - 10 7 Hz) arises from 
a Maxwell- Wagner relaxation caused by the membrane 
of the red blood cells [7l|. The 7-dispersion (at about 
10 10 Hz) is caused by the cooperative reorientational mo- 
tion of water molecules [3,|7l|. At v < 10 5 Hz, strong 
EP effects show up in the blood spectra of Fig. |51 The 
complete broadband spectra are well fitted with a single 
Cole-Cole DRC circuit (eq. ((H) with (3 = 1), connected 
in series to the bulk impedance. The latter is modeled by 
the sum of two "normal" Cole-Cole functions (i.e. they 
are defined in £*; eq. ([3]) with f3 = 1) to account for the 
f3- and 7-dispersions. The finding that a single DRC cir- 
cuit is sufficient to account for the data may be ascribed 
to the fact that the onset of the EP effects in blood oc- 
curs at rather low frequencies, comparable to the results 
in CKN (Fig. |2J). 

The diamonds in Fig. [B]show the results for a 0.009% 
(5 mmol/1) Lysozyme solution. Protein solutions usu- 



ally contain free ions and indeed EP effects show up 
in the present spectra, leading to strong dispersion in 
both quantities. The very gradual transition to a low- 
frequency plateau in e'(v) at v < 10 2 Hz indicates the 
presence of a second, smaller step on top of the main dis- 
persion. Indeed, excellent fits of the experimental spec- 
tra were achieved by using two DRC elements in series to 
the bulk. The DRC circuit leading to the smaller, low- 
frequency step is Debye-like (a = 0, j8 = 1); the broad 
main dispersion step is characterized by a — 0.46 and 
(3 = 1. 

Finally, Fig. [6] also provides some typical data on elec- 
tronically conducting materials, where the EP arises from 
the formation of Schottky diodes at the electrode/sample 
interface. Even when the polarity of the oscillating field 
changes, one of the two oppositely poled diodes at both 
sides of the sample always is blocked and dominates 
the response at low frequencies. In recent years, es- 
pecially EP effects in non-metallic transition-metal ox- 
ides have gained increasing significance: Many of these 
materials show effects like magnetocapacitance, multi- 
ferroicity or the occurrence of very large dielectric con- 
stants that have triggered numerous dielectric investi- 
gations [HI, [72], [73|. However, EP effects can consid- 
erably hamper the detection of their intrinsic dielectric 
properties and even be misinterpreted as being of intrin- 
sic origin [24|, [2^, [zH- The examples provided in 
Fig. |6]include single-crystalline LagCuCvi, a parent com- 
pound of high-T c superconductors 21], single crys- 
talline Lai 5 / 8 Sr 1 / 8 Ni04, a charge-ordered CDC material 
(53j . and polycrystalline Pr 2 /3Cu3Ti40i2 [3l|, a material 
structurally closely related to the well-known CDC ma- 
terial CaCu3Ti40i2 [7(| . For all these materials, strong 
EP effects are revealed in Fig. [51 which look qualita- 
tively similar to those in the ionic conductors. In all 
cases, the spectra can well be fitted by using a DRC cir- 
cuit to account for the EP (lines). For La 15 / 8 Sr 1 / 8 Ni04 a 
second DRC circuit was used to cover the second disper- 
sion step seen below about 10 6 Hz. In ref. [5J|, by mea- 
surements with different contact types the high-frequency 
step was clearly identified to arise from EP. The origin of 
the low-frequency step is not clarified yet but it was spec- 
ulated that it may be connected to the inhomogeneous 
charge distribution in this prototypical charge-ordered 
material (53|. We also found the DRC to be applica- 
ble for the description of EP effects in Cu2Ta40i2 [H[ 
and CaCu 3 Ti40i2, another CDC material [z3]. 



5. SUMMARY AND CONCLUSIONS 

In the present work, we have provided dielectric spec- 
tra showing strong EP effects measured in a variety of 
materials. They partly belong to very different material 
classes like salt solutions, ionic liquids and melts, biologi- 
cal materials, and electronic semiconductors. Our results 
are intended to provide a comprehensive data base for 
comparing current and future approaches for the mod- 



8 



eling of EP effects in dielectric spectra. For this pur- 
pose, the data are made available in electronic form in the 
supplementary information. In addition, in the present 
work we have demonstrated that the DRC equivalent cir- 
cuit is an easy-to-apply phenomcnological model for fit- 
ting EP effects in dielectric spectra, which works well 
for all the different material classes investigated. How- 
ever, this model certainly cannot be regarded as an alter- 
native to microscopic models as discussed, e.g., in refs. 

El El EHU- 0n the other hand > 1* has proven 
useful for a quick, phenomenological description of EP 
effects, enabling fits of experimental data in the com- 
plete investigated frequency range, which can help in the 
determination of the intrinsic bulk properties of the in- 
vestigated material. 

In all the investigated ionically conducting materials, 
we find evidence for a second charge-transport process 
that leads to a second plateau in the frequency-dependent 
conductivity at frequencies below the commonly observed 
decrease of o'(y) due to blocking. In addition, in most 
of these materials we find a further decrease of o' (v) at 
the lowest frequencies investigated, indicating final block- 
ing of the charge carriers involved in this second pro- 
cess. These seemingly quite universal features have been 
only rarely detected so far: In most works on ionic con- 
ductors only the onset of the EP effects is observed in 
the spectra, and the data can be fitted by simple mod- 
els like an undistributed RC circuit or a CPE. It is un- 
likely that the second dispersion, observed in measure- 
ments extending to sufficiently low frequencies, is due to 
a second, much slower ion species or due to electronic 
conduction. Instead, to us an explanation in line with 
the model proposed in ref. [HI seems more likely: The 
same ion species moves fast in the bulk region and much 
slower in the space-charge region, thus leading to the 
observed two dispersion effects. In ref. 15;], a specific 
hopping-conductivity model [78| was used to model the 



transport process in the region close to the electrode. 
Here we use the DRC circuit with (3 — 1, which is equiv- 
alent to a CPE connected in parallel to a resistor [5fj| . A 
CPE, sometimes termed "universal dielectric response" 
[55| . is often employed to describe hopping conductivity 
in various materials and indeed it leads to very similar 
frequency dependence as the model considered in [l5| . 
Thus, using the DRC circuit is in line with the model of 
ref. [15[ and adding a second DRC element also accounts 
for the final blocking of the slow ions in the BE region. 

It is an interesting finding that the DRC element also 
provides a perfect description of the EP effects observed 
in electronic conductors. In fact, La2Cu04 (Fig. [6]) was 
the first material, where we introduced this circuit for the 
description of an electronic conductor 15 years ago [2lj |. 
Great care has to be taken to avoid the misinterpretation 
of relaxation-like features in the spectra of semiconduct- 
ing materials as intrinsic effects that could be indicative, 
e.g., of ferroelectricity or "colossal" dielectric constants. 
Thus the correct modeling of these effects is essential. 
The fact that the present DRC-circuit approach leads 
to a different frequency dependence than an intrinsic re- 
laxation process, modeled, e.g., by a Cole-Cole function 
(see, e.g., Fig. [3]), provides some possibility to distin- 
guish between intrinsic and extrinsic effects. In addition, 
of course the procedures proposed in ref. [24] (e.g., per- 
forming measurements with different contact types, etc.) 
should be performed to arrive at a definite conclusion. 
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